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I. Bausmerth 1 A. Rccati 1,2 , and S. Stringari 1 
1 Dipartimento di Fisica, Universita di Trento and CNR-INFM BEC, 1-38050 Povo, Italy 
2 Physik Department, Technische Universitat Miinchen, D-85747 Garching, Germany 

Using many-body results available from diagrammatic and ab initio Monte Carlo calculations we 
analyze the phase diagram /i = (/if + /xj.)/2 versus h — (fi-j — /ij.)/2 of a unitary Fermi gas at zero 
temperature with population imbalance and unequal masses. We identify different regions where the 
gas is superfluid, partially polarized or fully polarized and determine the corresponding coexistence 
conditions. The asymmetry in the phase diagram, caused by the mass imbalance, and its effect on 
the Chandrasekhar-Clogston limit for the critical polarization are explicitly discussed. The equation 
of state of the superfluid and normal phases is employed, within the local density approximation, 
to predict phase separated configurations in the presence of harmonic trapping potentials. 
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I. INTRODUCTION 

The study of polarization effects in Fermi superfluids 
has been the object of intense experimental and theo- 
retical work (for recent reviews on the subject see, e.g., 
[l], 0]) in ultracold atomic gases in the last few years. Cru- 
cial goals of these studies are the identification of quan- 
tum phases and the determination of the Chandrasekhar- 
Clogston limit of critical polarization above which the 
system is no longer superfluid. 

The recent observation of hetcronuclear Feshbach res- 
onances in ultracold mixtures of two fermionic species 
H as well as the realization of a degenerate two-species 
Fermi-Fermi mixture Q has opened new stimulating per- 
spectives in the field of Fermi superfluids built with 
atomic species of different masses. 

The phase diagram of Fermi mixtures with unequal 
masses and the corresponding polarization effects, includ- 
ing the possible occurrence of the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase Q, have already been the 
object of theoretical predictions based on BCS mean- field 
theory @, 1, 1 G3 . This theory is known to give rea- 
sonable predictions at unitarity in the case of unpolarizcd 
configurations (sec, for example, 0)- However, it fails to 
give quantitatively correct results in the imbalanced case 
and to predict the Chandrasekhar-Clogston limit of crit- 
ical polarization. The failure of the BCS mean-field is 
mainly due to the fact that it ignores the role of interac- 
tions in the normal phase which are now understood to 
play a crucial role at unitarity [ll|, E3] ■ 

The main goal of this paper is to use the present knowl- 
edge of the equation of state of Fermi mixtures with 
unequal masses to give quantitative predictions for the 
phase separation between the normal and superfluid com- 
ponents. Our analysis is based on the study of the zero 
temperature fj,-h phase diagram as shown in Fig. [T] of 
the uniform two component gas, where fj, = (/if + /if)/2 



is the chemical potential and h = (/tf — /tj.)/2 is an ef- 
fective magnetic field. The phase diagram at unitarity 
is determined thanks to the knowledge of the equation 
of state available from diagrammatic techniques applied 
to highly polarized configurations and from Monte Carlo 
simulations. The phase diagram is then used, in the local 
density approximation (LDA), to calculate the density 
profiles of the two Fermi components in the presence of 
harmonic trapping. 
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FIG. f: (Color online) In the case of equal masses the /t- 
h phase diagram is symmetric with respect to zero effective 
magnetiv field h. Shown are the superfluid (S, red), the par- 
tially polarized (PP, green) and fully polarized (FP, blue) 
phases. 



We begin in Sec. |TT] by reviewing the general theory 
of the normal state for equal masses, and discuss its ex- 
tension to the unequal mass case. Then in Sec. 11111 we 
introduce the bulk phase diagram and discuss its prop- 
erties in dependence on the mass ratio k. In Sec. IIVI 
we use the phase diagram in local density approximation 
to calculate the density profiles, focusing on three par- 
ticular configurations. Finally in Sec. |V] we draw our 
conclusions. 
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II. NORMAL STATE OF A FERMI GAS WITH 
UNEQUAL MASSES 

The equation of state of the normal phase in the uni- 
tary limit of infinite scattering length and at zero tem- 
perature was first derived by Lobo et al. [ll[ m the case 
of equal masses using the concept of quasiparticles. As a 
function of the concentration x = nj,/ n T and for i<l 
it is given by 
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where N-\ is the total number of spin-f atoms and Ep^ = 
?i 2 /2m(67r 2 7i|) 2 / 3 the Fermi energy of the spin-| gas. In 
Eq.([T]) it is assumed that adding a few spin- J. particles 
to a nonintcracting spin-f sea, the |-atoms form a Fermi 
gas of quasiparticles with an effective mass m* . 

The first term in Eq. ([1]) corresponds to the energy per 
particle of the noninteracting gas, while the term linear 
in x gives the binding energy of the spin- J, particles to the 
spin-f sea. The interaction between f and J, particles is 
accounted for by the parameter A = 0.99(1) [13j . which 
is proportional to the ratio The Fermi gas of 

quasiparticles with an effective mass m* contributes to 
the total energy in Eq. (TTJ) by the quantum pressure term 
proportional to cc 5 / 3 , and in the case of equal masses 
m*/m = 1.09(2) [HI]. Eventually, the term proportional 
to x 2 can be interpreted as an interaction between quasi- 
particles, and its value B = 0.14 has been determined fit- 
ting the expression (p} to the Monte Carlo results for the 
equation of state as a function of the concentration [l3T ]. 
Although based on a small x expansion, Eq. ([T]) turns out 
to account for the x-depcndcncc of the equation of state 
also for values of x close to 1. 

For m-f 7^ mj the values of the parameters A and m* 
depend on the mass ratio in a non-trivial way and have 
been calculated in [lj[ and [l~5| as functions of the mass 
ratio m|/?rij = k using diagrammatic many-body tech- 
niques. The parameter A is an increasing function of the 
mass ratio k going to infinity for k — > and reaching 
the asymptotic value A ~ 0.45 for k — > oo. On the other 
hand at unitarity the effective mass, which we will denote 
in the rest of the paper as m* = mj_-F(re), shows a weak 
dependence on the mass ratio. 

The quasiparticle interaction B has up to now only 
been determined for equal masses m j = m-f . We can find 
a first estimation for B(k) in the following way. In the 
unpolarizcd case (x = 1) the energy of the normal state 
as a function of the mass ratio k has been calculated using 
Monte Carlo methods 

nana, 

resulting in the expression 
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where the dimensionless parameter £n(k) accounts for 
the interactions, and m K = (m^rrii) / (m^ + mi) is the 
reduced mass. First results based on quantum Monte 
Carlo (QMC) calculations suggest that the dependence 



of the interaction parameter £n(«0 on the mass ratio k 
is very weak [l8| so that we can set £n( k ) = Cn = 0.56 

HE 

Il7j | also for k ^ 1. Therefore the effect of unequal 
masses on the energy of the unpolarized normal state en- 
ters only through the reduced mass m K . We can define 
B as a function of the mass ratio by requiring that the 
energy of the normal state be reproduced by Eq.|T]) for 
a concentration x = 1 with the K-dependent interaction 
parameters A and m* /mi given in hj. Then the gener- 
alization of Eq. |T]) to the unequal mass case yields 
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III. PHASE DIAGRAM 



In terms of the mass ratio k the superfluid energy takes 
the form 
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.(67T 2 n s ) 2 / 3 = e s (n s , K ), (4) 



where Ng is the number of atoms in the superfluid phase, 
ns the superfluid density, m K the reduced mass, and 
£s(k) accounts for the interactions in the superfluid. Also 
in the superfluid phase the coefficient £s(£L) has only a 
very weak dependence on the mass ratio [18[ so that we 
can set £s( K ) = £s = 0-42 [H, [ItJ as in the equal mass 
case. 

In order to establish the phase diagram for the system 
we address the equilibrium conditions for the phase sepa- 
ration of the superfluid and normal state in the bulk. We 
start by writing down the energy of the system at zero 
temperature 



/ dr 

dr 



e s (n s , n)n s - ^n s 
en(x,K)n<t - fx°n^ > ( 5 ) 



where es(«Sj K ) and e^(x, k) are the energy densities per 
particle, ng, and the densities, ^ and fi® the chemi- 
cal potentials of the spin-f and spin- J, component, respec- 
tively, and /ig = + /i < j')/2 is the superfluid chemical 
potential. 

To find the equilibrium conditions we minimize the en- 
ergy with respect to the densities of the superfluid and 
normal phase, and we find the chemical potentials 
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where prime means the derivative with respect to x. 
Eventually requiring that the pressure of the two phases 
be the same yields 
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Making use of Eqs.Q and ([3]) we can write the equal 



(9) 




0.1 0.2 



FIG. 2: (Color online) Critical concentration x c (k) for the 
bulk system to phase separate as a function of the mass ra- 
tio /t (solid blue). In comparison, also the concentration de- 
rived from the BCS mean-field solutions at unitarity is shown 
(dashed red). 

pressure condition as the density jump 



From Eqs. 
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we obtain the implicit equation 
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which gives the Chandrasekhar-Clogston limit x c in de- 
pendence on the mass ratio n. For values smaller than 
x c (k) the system remains normal, while for x > x c (k) the 
system starts nucleating the superfluid and phase sepa- 
rates into those two states. In Fig. [2] we plot x c (k) for 
mass ratios 0.1 < K < 10 (blue solid line). Compar- 
ing with x c (k = 1) = 0.42 we find that for mass ratios 
K > 1 the concentration needed to create a superfluid 
phase decreases, while for mass ratios k < 1 it first in- 
creases and reaches a maximum value at k ~ 0.3. In 
the same figure we plot x c (k) as resulting from the BCS 
mean-field approach at unitarity (red dashed line, see 
also e.g. 0). In the latter treatment interactions in 
the normal phase are not taken into account, and hence 
its energy is just the sum of the f and J, components, 



E® cs = Ep^N^ + Ep^Ni, and the interaction parame- 
ter for the superfluid is £f cs = 0.59. The significant 
quantitative difference between the two curves proves the 
importance of interactions fl2| . 

It is worth noticing that x c is sensitive to the actual 
value of the parameters used in Eq.Q. Since an exact 
calculation of the parameter B(n) in the case of unequal 
masses is still lacking, the interpolated value of B(k) 
might be a significant source of error. As for k >> 1 
the kinetic energy of the quasiparticles becomes irrele- 
vant (see Eq.Q), the Chandrasekhar-Clogston limit is 
only determined by the values of A(k) and B{k). Thus 
an uncertainty in B affects more our predictions. On the 
other hand for k > > 1 a polarized superfluid phase might 
have to be included in the description of the system so 
that the two-phase assumption will not longer be valid 
(see also discussion in Sec. IV) . 

We varied the value of B(k) by ±10% to see its final 
impact on the value of x c (n), and we find that the varia- 
tion in x c (k) is around =f5% for mass ratios k < 1, while 
for k > 1 it is F-10%. 

In terms of the chemical potentials of the | and j com- 
ponents the phase transition is characterized by the crit- 
ical value T) c (k) = (fJ-i/ fJ-^)x c - From the knowledge of 
f7c(ft) we are able to determine the coexistence lines be- 
tween the superfluid and the normal phase. 

We represent the different homogeneous phases em- 
ploying the \x — h phase diagram, where 2/i = ^ + fj,^ 
and 2h = /ij — fi^. The transition line between the su- 
perfluid (S) and partially polarized (PP) phase is given 
by 
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and stands for the first-order phase transition between 
the unpolarized superfluid and the partially polarized 
normal phase. 

The second-order phase transition between the par- 
tially polarized (PP) and the fully polarized (FP) phase 
occurs at x = 0, which corresponds to ^i/fJ.-\ = 
—3/5A(k), and thus the coexistence line is given by 
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Finally, the transition line between the fully polarized 
gas and the vacuum is given by the simple K-independent 
relation 



MFP 



(14) 



The phase diagram for unequal masses is not symmet- 
ric with respect to zero effective magnetic field h as can 
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be seen in Fig. [3l where we choose k = 2.2 correspond- 
ing to the case of a 87 Sr - 40 K mixture [l9||. While the 
supcrfluid S moves clockwise (anticlockwise) for k > 1 
(k < 1), the partially polarized PP moves in the oppo- 
site direction, see e.g. Figs. [T] and [3] In all the figures 
we use a solid line for the first-order phase transition, 
a dashed line for the second order phase transition, and 
a short-dashcd-long-dashed line for the transition to the 
vacuum. 

Such an asymmetry in phase diagrams is general for 
this system and has been already noticed by Parish et 
al. |9j in the T/llvs h/u phase diagram, and by Iskin 
and Sa de Melo in the P = (iV T - iV;)/(iV T + N t ) vs 
(l/kF,+,cLp) diagram. 
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FIG. 3: (Color online) For k = 2.2 the phase diagram is 
asymmetric. Shown are the superfluid S (solid red lines), 
partially polarized PP (dashed green) and fully polarized FP 
(dot-dashed blue) regions. 



In particular, from Eq. (|12p we can identify a critical 
mass ratio k* = 2.73 above which the superfluid region 
has moved entirely to the h > plane (see e.g. Fig. 
Hla)). This shift of the superfluid region above a certain 
mass ratio k has also been identified by Parish et al. @ 
applying BCS mean-field theory yielding n* ~ 3.95. 

At the same time, for k > k* the sum of the spin-f 
and spin- J, densities in the partially polarized phase is 
bigger than the superfluid density, (rij +rii) > 2ns- This 
anticipates the fact that in a trap the heavy partially 
polarized phase can sink towards the center, while the 
supcrfluid will form a spherical shell around it even if 
the two species feel the same trapping potential. This 
peculiar formation of a "sandwiched" supcrfluid has been 
previously identified also in @, H, 0] • 



IV. TRAPPED GAS 

Having constituents with different masses and hence 
different magnetic and optical properties permits to en- 
gineer different configurations in the trap depending on 
the mass ratio, the polarization, and the choice of the 
trap parameters. 



In order to study the trapped case we assume that the 
external potential is harmonic of the form V a {r) = ^a a r 2 
where a„ = m a ui 2 with a =T,|, and that the local 
density approximation is applicable. Thus the config- 
uration in the trap is found by using the expression 
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leading to 




(15) 



where we define 770 = mJ/a 4 ! as the central imbalance of 
the system, and (i?|) 2 = 2u^/a-\. Note that if aj = aj. 
the effective magnetic field h does not depend on the 
position in the trap but is only a function of the central 
imbalance 770. Concerning the central imbalance of the 
chemical potentials we have that if 770 < Vc{ K )i there is no 
superfluid and the system consists only of the partially 
and fully polarized component. In the case that ?yo > 
?7c(k), we have a superfluid component whose fraction is 
determined by the value of 770. 

In the following we will describe three different cases 
with different values of the polarization 



P 



(16) 



where the interplay between the asymmetry in the masses 
and in the trapping potential gives rise to configurations. 



A. Unequal masses with equal trapping 

We first analyze the situation when the spin-"f and 
spin- J. components have different masses k / 1 but feel 
the same restoring forces a-f = a±. This would be the 
case, for example, if the fermions are trapped magnet- 
ically and have identical magnetic moments. For equal 
populations and for mass ratios in the range 0.36 < K < 
2.73 the system is completely superfluid. In the opposite 
case, the system can never be completely supcrfluid even 
if the populations are equal. Therefore, we can also have 
the particular configuration of a system consisting only of 
a partially polarized phase (without any fully polarized 
part). 

If 770 > 1/ 7 7c(1/k) the trapped system will consist 
of a three-shell configuration, where the superfluid is 
sandwiched between a "heavy" normal phase (heavy spin- 
l arc the majority) at the center of the trap, and a 
"light" normal phase (light spin-f are the majority) in 
the outer trap region. 

As an example we choose the mass ratio k = 6.7 cor- 
responding to a 40 K - 6 Li mixture 0, [l^. The phase di- 
agram of the system is shown in Fig. |4ja) together with 
the LDA line for a central imbalance 770 = 0.3 (black 
dot). The intersection of the LDA line with the coexis- 
tence lines determines the radii of the configuration, from 
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FIG. 4: (Color online) (a) Phase diagram for n = 6.7, cor- 
responding to a 4() K- 6 Li mixture, and the LDA line (black 
vertical) for a central imbalance r/o — 0.3 (black dot) and 
«l = ax. (b) Density profiles for P = —0.13; the inset shows 
a zoom into the outer superfluid-" light" normal border. 

which we are able to calculate the density profiles. These 
are shown in Fig. BJb) for a polarization P = —0.13. 

The density jump (or drop) between the superfluid 
and both normal phases is a function of K according to 
Eq. (fT0"|) . For k = 6.7 at the "heavy" normal - superfluid 
border, n.j ~ 1.92ns and n-f = x c (—)nx ~ 0.86ns, while 
at the superfluid- "light" normal border nj ~ 0.71ns and 
nx = x c (n)nx ~ 0.17ns (see inset Fig. H{b)). Note that 
this is quite different compared to the equal mass case, 
where the jump between the superfluid and the majority 
component is n-\ ~ 1.01ns and hence the spin-f density 
is practically continuous. 



FIG. 5: (Color online) (a) Phase diagram for k = 2.2 and 
LDA line for 770 = 2.1 (black dot) and aja^ = 8. With 
this choice, it crosses the "heavy normal", superfluid, "light 
normal" and fully polarized phases, (b) Density profiles for 
a global polarization P — 0; the inset shows a zoom into the 
superfluid-" light" normal border. 

k c ~ 6.7. In Fig. EJa) we show the fi — h phase dia- 
gram of such a particular configuration, where we choose 
k = 2.2 corresponding to a 87 Sr - 40 K mixture [l9(. The 
LDA line is drawn for the values rjo = 2.1 (black dot) and 

"f/ a T = 8 - 

In the density profiles as shown in Fig. EJb) we have 
chosen the parameters such that the resulting global po- 
larization is P = 0. 



C. No trapping for | component 



B. Unequal masses with trapping anisotropy 

Using unequal restoring forces for the trapped atoms 
the mass ratio for having a sandwiched superfluid needs 
not to be necessarily bigger than the critical value k* . In 
order to have a three-shell configuration the condition is 
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K 7a > 
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(17) 



For example, for equal trapping frequencies, Eq. p7|) sim- 
plifies to k > [l/?7 c (l/re)] resulting in the critical mass ra- 
tio k c ~ 1.95, while for equal oscillator lengths one gets 



An interesting limiting case is when one of the clastic 
constants a a is zero (or very small), implying that one 
of the components would not be confined in absence of 
interspecies atomic forces. 

If we assume that aj — > 0, the LDA line in the fi-h 
phase diagram is parallel to the polarized-vacuum tran- 
sition line as shown in Fig. [^a). 

Let us start considering the equal mass, nif = mj = 
m, highly unbalanced JVj. 3> JVf case. The densities are 
easily found to be [l^l 



[^nx{v)f' 3 + Vx{v\ 
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FIG. 6: (Color online) (a) Phase diagram for k — 1. The 
LDA line represents the case Qf = 0. (b) Density profile in 
the limiting case a-\ = and P = —0.42. 

M?' = £,[6^(r)f 3 + V{(v), (18) 

where /if = ^ + §^\ V-[{r) = ^(rHfAV^r) andA = 
A(k = 1). From these equations it is clear that if Vj — * 0, 
the "[-atoms feel nevertheless the renormalizcd potential 
|j4V|(r) and are confined due to the interaction with the 
|-component. In this regime /i^ is negative and in the 
limit of a single f-atom, i.e. /i| — > 0, it takes the value 
/i| = — corresponding to a polarization P = — 1. 

This induced trapping mechanism would not be predicted 
by a BCS mean-field description, where interactions are 
absent in the normal phase, and the "["-atoms cannot be 
confined by the |-atoms. 

Increasing the number of "[ particles, the LDA line 
moves upward until it crosses the origin of the phase 
diagram, corresponding to (j^ = 0, and the system re- 
mains normal since for equal masses the slope of the 
superfluid-partially polarized coexistence line is bigger 
than the slope of the LDA line, i.e. t] c (1/k — 1) > 
(see Eq. ((T2|) ). Moreover, in this case the |-fully polar- 
ized phase is absent as the radii of the J. and "[ species 
coincide, and in this limit the polarization approaches 
the value P = -0.42. 

If we further increase we enter in a three-shell con- 
figuration including an intermediate superfluid compo- 
nent. But since in this case the atoms of species "[ are no 
longer confined, they escape from the trap, and the sys- 



tem goes back to the normal state previously described. 
Hence we can never find a stable configuration contain- 
ing a superfluid region, and the polarization of the system 
will always be in the range — 1 < P < —0.42. 

Note that the same scenario is valid for mj > 0.9to|, 
where T] c (1/k) is positive and the range of the polariza- 
tion is between P = — 1 and an upper value which is 
dependent on k. 

Interestingly, in the case mj, < 0.9mj, for which 
r) c (l/n) < 0, we find that adding "["-atoms we end up in 
a superfluid state (20| characterized by a density profile 
given by 

^=&^[67r 2 n s (r)] a/8 +y(r), (19) 

where V(r) = ^mio^r 2 is the effective potential felt by 
the superfluid. This configuration would correspond to 
a LDA line which stays entirely in the superfluid region, 
crossing the origin of the phase diagram. The value of 
the polarization for m[ < 0.9mf covers the entire range 
-1< P< 0. 

From the experimental point of view the above config- 
urations could be in principle reached starting with both 
the trapping frequencies different from zero and a certain 
initial polarization, and then opening adiabatically the 
trap for the "[-atoms. For instance, starting with only a 
superfluid in the trap the final state of the system will be 
simply a superfluid with a bigger radius for < 0.9m|, 
while it will be a normal state in which both components 
have the same radius (see Fig. [5{b)) if Ti| > 0.9m-f. 



V. CONCLUSIONS 

We have studied the zero temperature fj.-h phase di- 
agram of the unitary Fermi gas in the case of unequal 
masses, assuming phase separation between an unpo- 
larized superfluid and a polarized normal phase. The 
latter is described by an equation of state which, un- 
like in the BCS mean-field treatment, takes into account 
the effect of the strong interaction. As we have shown, 
this has a dramatic impact on the results such as the 
Chandrasekhar-Clogston limit needed to start nucleating 
a superfluid. 

Using LDA we have determined how the trapped con- 
figuration depends on the trapping potential, the mass 
ratio, and the polarization. Many different configura- 
tions are possible. Among them it is worth mentioning 
the three-shell configuration @, H, @, where the super- 
fluid is sandwiched between a "heavy" normal phase at 
the center and a "light" normal phase towards the edges 
of the trap. Note that the shells can occupy quite small 
regions, and we cannot exclude that surface tension plays 
an important role in this case. 

We can also have non-trivial configurations even if one 
of the two components is not trapped, but still remains 
confined due to the interaction induced trapping. Such 
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configurations can be experimentally obtained by adia- 
batically opening the trap for one of the two species. 

An important issue is the existence of other phases at 
unitarity. In the present work we assume that only two 
phases are possible, and hence we have not considered 
any polarized superfluid state. For the equal mass case 
the assumption seems to be correct and is theoretically 
understood by comparing the phase separated state en- 
ergy with the polarized superfluid energy calculated via 
Monte Carlo, as in e.g. [2l|. The same information is not 
yet available for the unequal mass case. However, taking 
the quasi-particle point of view in [2l| and the recent cal- 
culation for equal population by Baranov et al. [22j|, it 
seems that when the mass of the minority component is 
much bigger than the one of the majority component, the 
polarized superfluid phase should be included in the de- 
scription, as predicted by mean-field theory Theoretical 



work in this direction is in progress. 

Note added. Recently, a Monte Carlo calculation for 
the system considered here for k — 6.5 was posted in [23j |. 
Although an accurate comparison has not yet been done, 
it seems that the Monte Carlo analysis is in agreement 
with our description. 
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